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Abstract 

We study two-dimensional, large N field theoretic models (Gross-Neveu model, 
't Hooft model) at finite baryon density near the chiral limit. The same mechanism 
which leads to massless baryons in these models induces a breakdown of translational 
invariance at any finite density. In the chiral limit baryonic matter is characterized 
by a spatially varying chiral angle with a wave number depending only on the 
density. For small bare quark masses a sine-Gordon kink chain is obtained which 
may be regarded as simplest realization of the Skyrme crystal for nuclear matter. 
Characteristic differences between confining and non-confining models are pointed 
out. 



1) Introduction 



The description of baryonic matter on the basis of QCD remains a theoretical challenge, 
especially since lattice gauge calculations have so far been of little help for this problem. 
Nuclear physics, relativistic heavy-ion physics and astrophysics are some of the fields 
which would greatly benefit from any progress on this issue. Recently, a new surge of 
interest has been triggered by the suggestion that at high density, the novel phenomenon 
of color superconductivity might set in |[J §]. This development has highlighted how little 
is known reliably about dense, strongly interacting matter. 

Here, we address a much simpler finite density problem where a full analytic solution 
can be found: We consider two-dimensional model field theories with interacting fermions 
at or near the chiral limit. Specifically, we have in mind the two-dimensional version of 
the Nambu-Jona-Lasinio model ||, i.e., the chiral Gross- Neveu model Q, and QCD2 
with fundamental quarks, the 't Hooft model |5fl. In both cases, one considers a large 
number iV of fermion species and investigates the limit N —>■ 00, keeping Ng 2 fixed 
@. These two models are quite similar as far as their chiral properties are concerned 
but differ with respect to confinement of quarks which is only exhibited by the 't Hooft 
model. So far, the phase diagram of the Gross-Neveu model has been studied extensively 
as function of temperature and chemical potential |7j ||, ||, [10], [11], and the results seem 



to be uncontroversial. The 't Hooft model on the other hand has been investigated only 
sporadically at finite temperature ]20|, [2TJ, most recently in Ref. [^2[], but hardly anything 
is known yet about its properties at finite density |^T] . 



Our point of departure is the following observation: Both of these models possess 
light baryons whose mass vanishes in the chiral limit [23, |23|, |2J)| (by light, we mean that 
Mb/N is small on the relevant physical scale). This is of course no accident, but a generic 
feature of models with broken chiral symmetry in 1+1 dimension — the baryons are 
topologically non-trivial excitations of the Goldstone boson ("pion") field. By contrast, 
the Gross-Neveu model with discrete chiral symmetry can only accommodate baryons 
whose mass scales with the physical fermion mass and hence stays finite in the chiral 
limit. Nevertheless, it has been argued that both variants of the Gross-Neveu model have 



identical phase diagrams JTl|. Since we find it rather perplexing that the structure of 
the single baryon should have no influence on the structure of baryonic matter, we have 
reinvestigated this issue. We have found that a combination of large N techniques with 
the strong constraints arising from broken chiral symmetry is powerful enough to allow for 
a simple, analytic solution of this problem in the vicinity of the chiral limit. The results 
of our analysis differ qualitatively from the conventional wisdom about the Gross-Neveu 
model and carry over to the (confining) 't Hooft model as well. They seem to confirm a 
number of investigations of other field theories at finite fermion density, where strikingly 
similar behaviour was found. These include exact studies of two-dimensional models like 
the massive |T2[ and massless Schwinger model jnj or finite iV c and Nf massless QCD 2 
HI as well as more approximate treatments of 4-dimensional large N QCD [[15], [16], [lTj 



and effective chiral quark models |18l [L£ 



A word of caution is in order here: Throughout this paper, we shall constantly 
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deal with spontaneous breaking of continuous symmetries and Goldstone bosons in 1 + 1 
dimensions, in seeming conflict with the Coleman-Mermin- Wagner theorem P7[ . It 
is well understood by now that the large N limit enables one to circumvent this no- 



go theorem. As clearly explained by Witten |p8| , the bad infrared behaviour of the 
boson propagator, when exponentiated, gives rise to a power law correlator \x — y]^ 1 ^ 
which becomes constant in the limit N —>■ oo. This kind of almost long-range order 



is also familiar from the two dimensional Xl^-model [29, 30]. Alternatively, one may 
argue that the mean field approximation predicts symmetry breakdown and that this 
result is protected against fluctuations (which would otherwise restore the symmetry in 
two dimensions) by 1/N suppression factors. In this sense, low dimensional large N 
theories are not only more tractable, but also physically more appealing than their finite 
N counterparts. They bear more resemblance to the real, 3+1 dimensional world. 

This paper is organized as follows. In Sect. 2 we briefly review the conventional 
analytical treatment of the chiral Gross-Neveu model at finite density and point out a 
certain deficiency of this approach. In Sect. 3 we repeat a similar analysis for the 't Hooft 
model, supplementing the analytical methods by numerical computations where necessary. 
In Sect. 4 the Skyrme [31] type of approach to the light baryons in both models is recalled 
|24| and generalized to the case of baryonic matter in the strict chiral limit. In Sect. 5 



we then allow for a small symmetry breaking mass term and make contact with the sine- 
Gordon kink chain, the two dimensional analog of the Skyrme crystal |32| . This is followed 
by a short summary and conclusions in Sect. 6. 



2) Chiral Gross-Neveu model at finite density: Con- 
ventional approach 

Let us first recall the standard treatment of the Gross-Neveu model at finite density. In 
the large N limit mean-field techniques become exact. Technically, they may be phrased 
in a variety of ways. We choose the language of relativistic many-body theory, following 
Refs. f2~4] , |25| , [33[[ , which we find particularly intuitive for the problem at hand. Then 
the vacuum, the baryon and baryonic matter are all described by a relativistic Hartree- 
Fock approach (for baryons in the large N limit this was first recognized in Ref. [0]). 
"Conventional approach" in the title of this section refers to translational invariance - 
we shall assume that the system is described by an interacting Fermi gas with prescribed, 
homogeneous density. We shall first deal with the Gross-Neveu model with continuous 
chiral symmetry (-0 — > e iai5 ip) and Lagrangian density M 

£ = ^i# + ^ 2 ((^) 2 + (^75^) 2 ) • (1) 

As a matter of fact, the corresponding calculation would be identical for the model with 
discrete chiral symmetry only (ip — > 75VO, where the 75-term in Eq. (|l|) is omitted. The 
results presented here are well known, but our aim is to criticize them in a novel way. 

We denote the fermion density per color (or baryon density) by p# = Pf/n (p/: 
Fermi momentum). At the mean field level, the fermions acquire a physical mass m 
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which has to be determined self-consistently. The ground state energy density per color 
is given by 

£ n f A / 2 dk r- — m 2 

— = -2/ — Vm 2 + k 2 + (2) 

N J Pf 2tt 2Ng 2 V ; 

where A is an ultra-violet cutoff. The first term is just the sum over single particle energies 
for all occupied states (the Dirac sea plus all positive energy states with \p\ < pf), the 
second term the usual correction for double counting of interaction effects familiar from 
the Hartree-Fock approximation. To renormalize the theory, let us first consider the limit 
Pf — *■ 0, denoting the physical fermion mass in the vacuum by m , 

£ /-a/2 dk /—„ — — m, 

Minimizing £ with respect to tuq yields the relativistic Hartree-Fock equation 



-2/_ ^m 2 + k 2 + ^- 2 . (3) 



-o(l + ^ln^)=0. (4) 

Due to the similarity in structure between the relativistic Hartree-Fock approach and 
BCS theory |B5| , this is often referred to as the "gap equation" . The non-trivial solution 
(which has always lower vacuum energy) yields the relation 

Ng 2 A 

— — In — = 1 5 

7r m 

which teaches us how the bare coupling constant depends on the cutoff parameter, given 
mQ. Recall that the Gross- Neveu model shares with real QCD both asymptotic freedom 
and dimensional transmutation; these properties are contained in Eq. (|5]). Using this 
relation to renormalize the matter ground state energy density, Eq. (^), we find (dropping 
an irrelevant term — A 2 /8tt) 

i i ( ni -I- . I m 2 -X- n 2 \ 

N 4tt ' 2vr" J V-/ ' - o_ i ... / • !(,; 

The energy is minimal provided m satisfies 



1 m 



i.e., for 



/ 2pf ( mo\ 
m = or m = m \ I \Pf < ~^r) ■ (8) 



m V 2 



The corresponding energy densities are 
£ v) 



N 
£_ 

N 



m=0 



2tt ' 



£ ml p f m Q p) ( m 
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Pf<^ • (9) 
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The physical quark masses (|j) and the energy densities (|9[) are plotted in Figs. 1 and 
2. From these figures one might be tempted to conclude that chiral symmetry is broken 
at low densities and gets restored in a second order phase transition at pf = m /2. As 
is well known, this does not occur, rather there is a first order chiral phase transition at 
Pf = mo/ '\/2. This can easily be inferred by inspection of the thermodynamic potential of 
the Gross- Neveu model For our purpose, the following physical reasoning is perhaps 
more instructive: Let us compare the energy densities (0) with the energy density for 
a system of size L divided into two homogeneous regions I (size £) and II (size L — £). 
In region I chiral symmetry is restored; it contains the extra fermions needed to get 
the prescribed average density (the "MIT bag" @). Region II consists of the physical 
vacuum with broken chiral symmetry, void of excess fermions. The mean energy density 
obtained in this way is 

£ ( L- 



N 

Minimization with respect to i yields 




V2pfL 



(10) 



(11) 



m 

valid for pf < m /V2, and hence the optimal energy density 

£ _ ml p f m ( m \ 

As shown in Fig. 2, this solution is lower in energy than the homogeneous one; moreover, 
it yields the convex hull of £ . It ends exactly at the first order phase transition point 
Pf = m /v^2 where all space is filled with one big bag. This should be contrasted to 
the scenario underlying Fig. 1 where the fermion mass decreases continuously. We thus 
recover the generally accepted mixed phase interpretation of the Gross-Neveu model at 
finite density. Notice also that only the total size of regions I and II matters, not how 
they are subdivided; there could be baryon "droplets" as well. Alternatively, the mixed 
phase curve in Fig. 2 with its linear dependence on pf could have been inferred from a 
standard Maxwell construction. It is interesting that a very similar qualitative behaviour 
was found recently in 3+1 dimensions, where the close relationship with the bag model 
was also stressed 0. 

One important point to which we would like to draw the attention of the reader is the 
behaviour of £ near ps = pf/ir = 0. Since ultimately, at very low density, the fermionic 
matter problem must reduce to the problem of a single baryon, one would expect 



d£ 



dp 



B 



= M B (13) 

PB=0 



where Mb is the baryon mass. In the present calculation, Mb is not the physical baryon 
mass, but the mass of an alleged "delocalized" baryon. This is inherent in the translation- 
ally invariant Hartree-Fock approach, i.e., the assumption that the single particle orbitals 
are momentum eigenstates. Using Eq. (O) we obtain in the homogeneous, single phase 
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calculation, Eq. (^|), M B = Nm , consistent with a short range force and a delocalized 
baryon. The (physically more viable) mixed phase approach, Eq. (T2), predicts a baryon 



mass lower by a factor of l/y/2. This factor can readily be understood in terms of the 
bag model. Indeed, it follows from Eq. fllPf) that (E = £L) 

'•-*-«@ + *)' (14) 

For B — 1, this expression can be interpreted as energy of a single baryon, I being its 
diameter. The first term is just the bag pressure (the difference between the energy 
density of the physical vacuum and that of the perturbative one, cf. Eq. (H) for pf = 0), 
the second the kinetic energy of iV massless quarks. The bag size £ is found through 
minimization of the energy (for B — 1) to be 

t=^-. (15) 
m 



Inserting this result into Eq. (H), one finds that the bag pressure and the quark kinetic 
energy contributions are exactly equal in this model and that Mb = Nitlq / a/2. 

However, the Gross- Neveu model possesses bound baryons with lowest mass Nmo/7r 
(kink solution for the model with discrete chiral symmetry [|], Q), or even massless 
baryons (model with continuous chiral symmetry p4j|). These binding effects are not 



1/N suppressed and should be correctly reproduced in a Hartree-Fock approach, in the 
low density limit. They have obviously been missed here due to our tacit assumption of 
translational invariance. There is no good reason why such effects should not play a role 
at higher densities as well. Moreover, differences between the continuous and discrete 
chirally symmetric Gross-Neveu models based on their different baryon structure and 
masses are not at all captured by the "conventional" approach. Below, we shall present a 
cure for this disease. Before that however, let us first repeat the naive calculation for the 
't Hooft model, where the corresponding results are not yet available in the literature. 



3) 't Hooft model at finite density, assuming transla- 
tional invariance 

The 't Hooft model is defined as the large N limit of 1+1 dimensional SU(iV) gauge theory 
with quarks in the fundamental representation ||, 

£ = ^ip^ - UiF^ . (16) 

Since the light-cone approach originally used by 't Hooft to determine the meson spectrum 
seems to be less convenient for the vacuum, baryon and baryonic matter problems, we 



shall work in normal coordinates. This approach was pioneered by Bars and Green [37 
and further developed in Refs. |38], ^ [25| . Common to all of these works is the fact that 
the gluons are gauged away (axial gauge), leaving behind a theory of fermions interacting 
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via a linear Coulomb potential. We refer the reader to the detailed derivation of the 
Hartree-Fock approach in Refs. P4"| , and immediately proceed to the formulae which 
are relevant for our purpose. Let us first summarize the treatment of the vacuum. A 
central quantity is the single particle density matrix in momentum space, 



Pip) = g + 7°Po(p) ~ VpiO) + 7 5 Ps(p) • 



(17) 



Its precise definition in terms of the quark fields is 



The Slater determinant condition characteristic for Hartree-Fock, 

p\p) = Pip) -> PUP) + Plip) + PUP) = 1/4 
holds manifestly in the parametrization 



( Poip) 
Piip) 
V Ps(p) 



/ sin 9(p) cos if 
sin9(p) sin ip 
\ cos9(p) 



(19) 



(20) 



9(p) is the Bogoliubov angle, (p the global angle which locates the broken symmetry 
vacuum on the chiral circle (hence it has no p-dependence) . The Bogoliubov angles are 
the variational parameters; this terminology stems once more from BCS theory, which 
has the same formal structure as relativistic Hartree-Fock. We choose (p — 0, the value 
reached if one lets the bare quark mass approach zero starting from a finite value. Then, 



Pip) = v{p)v\p) 
with the (positive and negative energy) Hartree-Fock spinors 



uip) 



cos 9(p)l '2 \ 
sin 6 ip)/ '2 ) 



v(p) 



-sin0(p)/2 
cos 6 (p)/ 2 



The vacuum expectation value of the Hamiltonian density reads 



£ 

N 



dp n . , No 2 rdp f dp' cos(0(p) — 0(p')) 
- pcoa 9{p)-—J- 



2tt 



2tt 



ip ~ P') 



(21) 



(22) 



(23) 



where the first term is the kinetic energy, the second the Coulomb interaction of the 
quarks. Varying with respect to the Bogoliubov angles 9(p), the gap equation is obtained 
in the form 



p sin 9 ip) + 



Ng 2 rdp'sin(9ip) - 9{p')) 



2tt 







(24) 



ip — P') 2 

The integral has to be defined as principal value integral, cf. Refs. |J7|, B?| |2"5|] . We shall 
also need the expression for the quark condensate in the vacuum, 



-N 



TP sin%) 



(25) 
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and the quark single particle energies, 



U (p) = p cos 6( P ) + —j — {p _ pf)2 . (26) 



Although the gap equation (|24| ) for the 't Hooft model must be solved numerically, the 
value of the quark condensate (|25|) is known analytically, owing to an indirect determina- 
tion via sum rules and the 't Hooft equation for mesons [JJ{J; it is 

N (Nq 2 \ 1/2 



The single particle energies (^6|) are badly infrared divergent, a source of a long and 
ongoing debate in the literature fllPf . To exhibit the divergence, we follow Ref. p5 



isolate the divergent part of the integral and regularize it by using a finite box of length 
L, 

( x n ( ^ N 9 2 f <ip' cos{e(p)-e(j>'))-l Ng 2 L 
W (p) = pcos%) + — ) — + — . (28) 

This last constant diverges for L —>■ oo but seems to be essential to account for confinement 
in such an independent particle picture: The isolated quarks behave roughly as if they 
had infinite mass. If one simply throws the infinite constant away, as is often done, one 
gets an awkward sign change in uj{p) at some low momentum p, cf. Refs. [EM 123, EH], and 



runs into serious inconsistencies in finite temperature Hartree-Fock calculations |21], |22|| . 
Fortunately, the constant drops out of the calculation of color singlet mesons, as already 
noticed by 't Hooft (his IR cutoff parameter A is related to our box size L by A = 12/ttL, 
cf. Ref. [pjj): The infinite self-energy term is cancelled by an equally infinite piece in the 
Coulomb interaction. We also note in passing that if one employs a finite box as infrared 
regulator, one is unambiguously led to 't Hooft's treatment of the quark self-energies rather 
than to Wu's alternative regularization prescription HTI, HO]. Since the emergence of the 



constant Ng L/48 in the single particle energy (p8|), but not in the vacuum energy (|23|). 
is rather important for our discussion and somewhat hidden in Ref. ||25|| , we have included 
a simplified version of the arguments underlying Eqs. fl23|) and ( ^8|) in the appendix. 

After this review of the treatment of the vacuum, we are in a position to include a 
finite baryon density, assuming translational invariance. If pf denotes the Fermi momen- 
tum, we have to replace the density matrix (|21f) by 

p(p) = Q(pf - \p\)u(p)u\p) + v(p)v\p) 

= e( Pf -\ P \) + e(|p| - p f )v{p) v \p) (29) 

where we have used the completeness relation for the spinors in the second step. In the 
expression for the Hartree-Fock ground state energy density (p3|) , according to the second 
line of Eq. (p9|), we must exclude the region [— p/,p/] from the momentum integrations 
and pick up an additional term due to the change in the baryon density trp, 

£ f dp 
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Ng 2 fdp fdp' ^ cos(%) -%'))-! 

Q{\p\-Pf)Q{\p\-Pf) — ~r- 2 



2n J 2ix (p — p 

Ng 2 f dp f dp' 



/^/^-WW-ftli. (30) 



4 J 2n J 2n J [P ~ P'Y 

This yields at once the following finite density generalization of the gap equation, 

n , , Ng 2 fdp',, .. , sm(9(p) - 6(p')) ,. , , , . 

psine(p) + j-JLe(\p'\ - Pf ) \p_ pf) , }) = o , (\p\>Pf), (31) 

whereas the condensate now becomes 

(W) = -N J %(\p\- Pf )sm6(p) . (32) 



The gap equation (|31| ) can easily be solved numerically for various pf. The resulting 
condensate is shown in Fig. 3. We find that it decreases monotonically with increasing 
density, disappearing at a critical Fermi momentum 

^"°- U7 (^-) 1/2 • (33) 

This behaviour is strikingly similar to the corresponding result for the Gross-Neveu model 
depicted in Fig. 1, again suggesting some phase transition with restoration of chiral sym- 
metry at high density. Here we are not able to go on and discuss whether we are dealing 
with a first or second order phase transition. The reason lies in the following problem: If 
we compute the energy density (|30|) for the 't Hooft model we discover that subtraction 
of the value at pf = is not sufficient to give a finite result. Unlike in the Gross-Neveu 
model, the difference is still IR divergent. To be able to proceed, we enclose the system 
once more in a box of length L. We then find that the divergence is due to the last term 
in Eq. (|30|) (the one which does not involve the Bogoliubov angles) which now contributes 
the following double sum to the energy per color, 



E 

N 



d - = lt^ £ (34) 
div lovr peI n ^ (p _ n ^j n 



Here antiperiodic boundary conditions for fermions have been employed in the box regu- 
larization, and correspondingly the interval / is defined in the following way, 

I=[—rif,rif] for B = 2n/ + 1 odd , J=[— n/ — l,n/] for B = 2nf + 2 even . (35) 

The result (|34|) is even more alarming than the non-convex behaviour of £ in the Gross- 
Neveu model, Fig. 2, due to its L-dependence. Adding quarks to the vacuum causes the 
energy to increase by an infinite amount in the limit L — > oo. Evaluating the double 
sums in Eq. (|3~4D for low values of B, we obtain information on the origin of this divergent 
behaviour. For B = 1 (/ = [0,0]) in particular, the calculated baryon mass (to leading 
order in L) is 
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This is the same relation as Mb = Nm in the Gross-Neveu model except that the physical 
fermion mass is replaced by the infinite constant Ng 2 L/ '48 characteristic of confinement, 
cf . Eq. ( p8| ) . For larger values of B Eq. (|3~4| ) does not simply yield multiples of the baryon 
mass (^), but one finds "interaction effects" of the same order of magnitude as the mass. 
As far as iV-counting is concerned, this is still in agreement with Witten's analysis of 



baryons at large N f34 |. However, since these delocalized baryons are presumably not 



very physical in the 't Hooft model, we refrain from further discussing these effects. 

Summarizing, the problems encountered in the Gross-Neveu model with transla- 
tionally invariant baryonic matter again show up in the 't Hooft model, although in a 
much more severe form. The physics reason is clear: In the Gross-Neveu model the cost 
of distributing N fermions over the whole space is governed by their physical mass; in 
the 't Hooft model, due to confinement of quarks, the corresponding quark effective mass 
diverges with the volume. On the other hand, it is known that both models do possess 
massless, delocalized baryons in the chiral limit. Evidently, this has to be accounted for, 
and we conclude that the naive, translationally invariant Hartree-Fock approximation fails 
miserably in describing the properties of baryonic matter. 



4) Massless baryons and baryonic matter in the 
Skyrme picture 

The existence of massless baryons in the chiral limit of the 't Hooft model has been 
demonstrated by bosonization [f23|, fHJ, variational pi] , and light-cone JZIJ techniques. 
These exotic objects are characteristic for 1+1 dimensional models with broken chiral 
symmetry and, as such, also present in the chiral Gross-Neveu model. A particularly 
illuminating derivation is due to Salcedo et al. pl| . These authors point out that the 
potential energy in such models is invariant under local chiral transformations, unlike the 
kinetic term which is only invariant under global ones. This led them to the following 
variational ansatz for the one-body density matrix of the baryon, 

p{x, y) = e lx{xh5 Pv {x - y)e- ix(y)75 . (37) 

Here p v (x — y) is the vacuum density matrix. If the vacuum breaks chiral symmetry, one 
can generate with expression ([H]) a new (exact or approximate) Hartree-Fock solution 
which breaks translational invariance but can carry non-zero baryon number. As shown 
in |24[], the baryon density is given by 

p B (x) = tr(p(x, x) - p v (0)) = -d x x{x) , (38) 

7T 

so that the baryon number coincides with the winding number of the chiral phase x( x )i 

B= / L dxl^= X(L) " X(0) eZ. (39) 

JO TV IX 

(Notice that x{L) ~ x(0) must be an integer multiple of 7r since otherwise bilinear fermion 
observables would no longer be periodic.) For this topological reasoning it is again recom- 
mendable to work in a finite box of size L. The topological interpretation of the baryon 
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number also agrees with exact results of Ref. [f2I| which were not restricted to the large 
N limit. In the absence of an explicit quark mass term, the ground state energy obtained 
from Eq. (|37| ) is 

E[p] = E[p v ] + N f L dx^(d xX ) 2 ■ (40) 

JO 2,71 

This result holds independently of the specific model, since the potential energy does not 
contribute to E[p] — E[p v ]. Differences between various models are of course still present 
in the vacuum density matrix p v in Eq. (|37|) but do not manifest themselves in the baryon 
energy. Minimizing E[p] with respect to \ yields the free (static) bosonic equation 

d 2 xX (x) = 0, x(L)=x(0)+*B (41) 

with the solution 

X (x) = ttB . (42) 

Here xq is a parameter which reflects the breakdown of translational invariance. The 
baryon density is x-independent (ps(x) = B/L) as follows more generally from axial 
current conservation in the chiral limit However, the scalar and pseudoscalar con- 
densates acquire a non-trivial x-dependence, 

{ipip) = (tpip) v cos (2ixB(x - Xq)/L) , 
(fas-ip) = -W)v sin {2irB{x -x )/L) . (43) 

Since fluctuations of x( x ) describe the massless Goldstone boson field, the baryon picture 
emerging here is very similar in spirit to the Skyrme model |3T|. The fact that the baryon 



is a topological soliton will become somewhat clearer once we include a small bare quark 
mass (see Sect. 5) but this solitonic character also holds in the strict chiral limit considered 
here. 

We can now discuss the baryon as well as baryonic matter from this point of view. 
The single baryon (B = 1) is spread out over the whole space, the chiral phase xi x ) 
making one turn with constant speed to minimize the kinetic energy (Fig. 4). The baryon 
energy is, using Eqs. (|40|442D , 



Eb = N^. (44) 

This confirms that indeed the baryon becomes massless in the limit L — > oo. Incidentally, 
expression (|44] ) is identical to the kinetic energy of N non-interacting, massless quarks in 
the lowest momentum state available for antiperiodic boundary conditions. Nevertheless, 
we are not dealing with the free, chirally symmetric theory, but with the broken phase of 
an interacting theory where the quarks are massive or even confined. 

It may be worthwhile to contemplate the structure of the baryon for a moment 
from the point of view of the relativistic Hartree-Fock approximation. In chirally non- 
invariant models one would suspect that the baryon comprises the filled Dirac sea plus 
one filled, positive energy valence level. This is exactly what one finds analytically in the 
non-chiral Gross-Neveu model |[33|| , or numerically in QCD 2 with heavy quarks p4|. The 
picture implied by the ansatz (|37|) in the chiral limit is rather different though. Denoting 
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the negative energy single particle orbitals in the Dirac sea by (p^"\x) (solutions of the 
first quantized Dirac equation with Hartree-Fock potential), the Skyrme type baryon ( p?TD 
admits the density matrix 

p(x, y) = e™^ L Y, ^l~ ) (x)^i- )t (y)e- i ^/ L (45) 

k 

where 

<pi~\x)=e ikx v(k) , (46) 

and the k are discrete momenta appropriate to the interval of length L. We first observe 
that the chiral phase factor splits the momenta of the right- and left-handed components 
into k±7r/L. Since the transformed single particle wave functions are no longer momentum 
eigenstates, translational invariance is lost. Secondly, we note that the presence of an 
infinite Dirac sea is crucial for getting the extra baryon charge, rather than a single 



valence state. If the sum over occupied states k in Eq. (45) was finite, we would trivially 
conclude that p(x, y) and p v (x — y) belong to the same baryon density (pb(x) = trp(x, x)). 
Due to the infinite number of occupied states however, p v (x — y) develops a singularity at 
x = y, and one has to do a more careful point splitting in order to compute the baryon 
density. The divergence is due to the UV region and therefore determined by the free 



theory (for more details, cf. Ref. [24]), 

lim tr(p(x, y) - p v (x - y)) = limtr |e^/ L (±5{z) - ^) - = \ . (47) 

The result 1/L is the baryon density for B = 1. This mechanism is strongly reminiscent of 
the calculation of anomalous current commutators, for instance in the Schwinger model. 
The extra baryon number does not reside in a valence level added on top of the Dirac 
sea but somehow emerges from the bottom of the Dirac sea if one modifies all the levels 
slightly — it is a vacuum polarization effect. 

Equipped with this exotic kind of baryon, we can now easily find the ground state 
of the system for any baryon density. As discussed above and illustrated in Fig. 4, the 
single baryon consists of one turn of a "chiral spiral" (parametrized by x( x )) over the 
total spatial length L of the system - - admittedly a somewhat elusive object in the 
thermodynamic limit. A finite density ps = B/L = Pf/n on the other hand implies that 

X(x) = p f (x - x ) , (48) 

i.e., one full rotation over a physical distance which has a well defined limit for L — > oo, 
namely The baryon density remains constant in space, but the condensates are 

modulated as 

(^) = (tpip) v cos 2p f (x - x ) , 
(V>i75?/>) = -{ijjijj) v sm2pf(x - x ) . (49) 

They can be viewed as projections of a "chiral spiral" of radius | (4>?p) v \ onto two orthogonal 
planes, see Fig. 5. This state breaks translational symmetry; it is a crystal. In fact, it 
may be viewed as the simplest possible realization of the old idea of a Skyrme crystal ||32||, 
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here in the context of large N two-dimensional field theories. One cannot tell where one 
baryon begins and ends — each full turn of the spiral contains baryon number 1. Only 
the condensates reveal that translational symmetry has been broken down to a discrete 
subgroup. The energy density of this unusual kind of "nuclear matter" is simply (after 
subtracting the vacuum energy density) 

Surprisingly, this is exactly what one would expect for a free Fermi gas of massless quarks 
although Eq. (p0|) holds for interacting theories where the vacuum has lower energy due 
to chiral symmetry breaking. In Fig. 6 we compare the energy density for this state to the 
ones discussed above for the Gross- Neveu model, where translational symmetry had been 
assumed. The crystal is always energetically favored, the dependence on pf is now convex, 
and there is no trace of a phase transition, neither first nor second order, at any density. 
The horizontal slope at pf = correctly signals the presence of massless baryons and 
eliminates the above mentioned problems with the spurious massive, delocalized baryons. 
We cannot even draw the corresponding picture for the 't Hooft model, simply because 
in this case the quark Fermi gas is infinitely higher in energy than the Skyrme crystal for 
L — > oo. Nevertheless, all the results for baryonic matter discussed in this section apply 
to the 't Hooft model as well. 

In the high density limit the oscillations of the condensates become more and more 
rapid. If we are interested only in length scales large as compared to 1/pf, the condensates 
average to zero. In this sense, one might argue that chiral symmetry gets restored at high 
density, although not in the naive way suggested by Fig. 2. 

Finally, we remark that the "chiral spiral" ground state for fixed baryon density still 
preserves one continuous, unbroken symmetry, namely the combination of translation and 
chiral rotation generated by P + PfQ§ (P: momentum operator, Q 5 : axial charge). One 
would therefore predict that RPA excitations on this ground state |24], |25| (or mesons in 
nuclear matter) will have only one collective, gapless mode, a hybrid of a "phonon" and 
a "pion". 



5) Non- vanishing bare quark masses 

In Ref. JMJ the Skyrme picture of the baryons in the 't Hooft model and chiral Gross- Neveu 
model was developed for small, finite bare quark masses, using the expression in Eq. ( |37D 
as a variational ansatz. For a single baryon, these authors have tested the accuracy of their 
procedure against the full, numerical Hartree-Fock calculation on a lattice. The results 
agreed perfectly at m q = 0.05 and were still surprisingly good at m q = 0.20, in units of 

\J Ng 2 /2-k. This makes it very tempting to speculate that the corresponding variational 
calculation can also give us a reliable picture of baryonic matter at finite density, away 
from the chiral limit. As compared to the formulae in the preceding section, the only 
change is the fact that the bare mass term now also contributes to the energy functional 
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Eq. 

E[p]=E[p v ] + N£dxl [ ^(d xX ) 2 + ^^ . (51) 

Here the condensate (^>ip) v refers to the vacuum in the chiral limit. Varying with respect 
to x( x ) then gives the static sine-Gordon equation f4*2]1 , 



2 2irm q {ipi(}) v . 
d x x + ^ sm 2 X = , (52) 

from which one reads off the "pion" mass (two dimensional version of the Gell Mann- 
Oakes-Renner relation El) 



2 A W)v , K a\ 

m n = -47rm ? — ^— . (53) 



The B = 1 baryon can be identified with the familiar kink solution of the sine-Gordon 
equation, 

X(x) = 2arctan(W a; - Xo) ) , (54) 



with mass 



M B = . (55) 

7T 



Since the single baryon has been discussed in detail in Ref. f24jl , let us immediately turn 



to multi-kink solutions as candidates for baryonic matter. Luckily, the sine-Gordon kink 
crystal has already been studied thoroughly in the literature, first in solid state physics 
t5[ and more recently as a toy model for the Skyrme crystal p6| , in terms of Jacobi 



elliptic functions and elliptic integrals ||47|| . We take over the results from Ref. |^6[ which 



is close in spirit to the present study although the authors did not have in mind two- 
dimensional large N field theories. Adapting the formulae of this work to our notation, 
the following steps allow us to generalize the Skyrme crystal of the previous section to 
small, finite bare quark masses: Let denote the mass of the Goldstone boson, Eq. (|53|), 
and pB — Pf/n the average baryon density (this is our definition of pf for the case of broken 
translational symmetry). We then first have to solve the transcendental equation 

= 2kK(k) (56) 

Pf 

for k where K(/c) is the complete elliptic integral of the first kind. The sine-Gordon kink 
crystal is then given by the following solution of Eq. (|52|) , 

7T Tfl 

X{x) = - + am(f, k) , £ = — -(a? - x ) , (57) 

(am(£, k) is the Jacobian elliptic amplitude function). From this, we can express the 
baryon density and the various condensates in terms of further Jacobian elliptic functions 
(dn, sn, cn) as follows, 

I TTl 

Pb{x) = -d x x{x) = — fdn(f, k) , 

7T Tlk 

W) = (^)vCos2x(x) = -(V^)v (cn 2 (£, k) - sn 2 (£, kfj , 
(-0175-0) = -(^) v sin2x(x) = (ipip) v 2sn(£, fc)cn(£, k) . (58) 
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Here £ is as defined in Eq. (|57|). Finally, the energy divided by the volume of this kind of 
matter is given by 

E(fc) denoting the complete elliptic integral of the second kind. 

Let us now illustrate these results in two regimes of interest, namely at low and 
high density. At low density (pj <C m^) k in Eq. ([56]) approaches 1 exponentially, and 
the baryon density features a chain of well resolved lumps whose shape is determined by 
the single kink solution (Fig. 7). Likewise, the condensates behave like those of a single 
baryon: (ipip) changes from the vacuum value outside the baryons to its negative in their 
center whereas (^i^ip) is peaked in the surface region of each baryon (Figs. 8, 9). These 
condensates are projections of the distorted "chiral spiral" shown in Fig. 10. The energy 
(|59| ) for low densities behaves as 

S « N^fl = M BPB , (60) 

showing the expected connection to the baryon mass. At high densities (pf ^> m n ), k 
approaches like 

ife«^. (61) 
Pf 

Thus £ in Eq. (|57|) becomes Pf(x — x ). Moreover, for k — > 0, the Jacobian elliptic func- 
tions am((, k), sn(£, k), cn(£, k) are known to reduce to the argument £ and the ordinary 
trigonometric functions sin£ and cos£, respectively. We thus recover the results for the 
simple chiral spiral in Sect. 4 (the parameter xq has to be readjusted to take care of the 
shift by 7r/2 in Eq. (0)). The energy in this case is approximately 

+ (62) 
iV 2rr 8n K ' 

The condensates look very much like the sin- and cos-functions of the massless case and 
need not be plotted. The baryon density wiggles around a constant value, reflecting the 
strong overlap of the baryons, and can be approximated at high density by 

Pb{x) w ^ U - 2~ (~^) sin 2 p/(x - x )j . (63) 

The behaviour of the baryon density Pb{x) as one increases pf (i.e., the mean density) is 
illustrated in Fig. 11. In the chiral- or high-density limit (m n /pf — > 0) Pb(x) eventually 
becomes ^-independent. This provides us with another way of understanding the struc- 
ture of matter described in the previous section, namely as arising from a chain of very 
extended, strongly overlapping lumps. 

Finally, let us come back to the question of validity of the variational calculation 
based on the chirally modulated vacuum density matrix (p7|), which we have left open so 
far. At very low densities when the individual baryons are far apart, we can presumably 
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rely on the numerical results of Ref. [p4 since the interaction effects between the baryons 



become small (as discussed in Refs. |24J and ^Bj , the baryon-baryon interaction is repulsive 
and falls off exponentially with the pion Compton wavelength) . At high densities, on the 
other hand, one would expect that a finite quark mass cannot make much difference as 
long as pf ^> m n . This takes us back to the massless case discussed in Sect. 4. In this 
limit, in turn, it is easy to convince oneself that the calculation becomes exact in the 
sense that one gets a true solution of the Hartree-Fock equation. Thus for instance for 
the 't Hooft model, the massless Hartree-Fock equation reads fl24 



u n (p^\x) = -i(-f 5 ) aP ^(pf\x) + J dy\x - y\p a p{x, y)<pf\y) ■ (64) 
Upon substituting 

^\x) = (e-^) a ^f\x) (65) 

as we are instructed to do by the ansatz P7|), we discover that does indeed solve 
the Hartree-Fock equation, the only change being that the single particle energy uo n gets 
replaced by u n +Pf. The same argument goes through in the chiral Gross-Neveu model, 
or in any field theory where the interaction term has a local chiral invariance. This proves 
that the result becomes exact in the chiral limit (to leading order in the 1/N expansion, 
of course) and makes plausible the hypothesis that it also correctly describes the high 
density regime for finite quark masses as long as pf 3> m w . 



6) Summary and conclusions 



In this paper, we have addressed the problem of baryonic matter in a certain class of ex- 
actly soluble field theoretic models, namely chirally invariant, large N, interacting fermion 
theories. We started out from a seemingly innocuous and well understood problem, the 
chiral Gross-Neveu model at finite density and identified one remaining weak spot: The 
energy density for baryonic matter in the standard Hartree-Fock approach does not have 
the correct low density limit which can be predicted from the known baryon spectrum 
of the theory. The origin of this problem which is not cured by a mixed phase approach 
conceptionally related to the bag model, is evidently the assumption of translational in- 
variance. Whereas this inconsistency can perhaps be ignored in the Gross-Neveu model 
(as it has been so far, to the best of our knowledge), in QCD 2 , it becomes fatal: Due to 
confinement, the analogous calculation yields an infinite energy for delocalized baryons 
or quark matter. This is unavoidable if one is careful in treating the infrared behaviour 
of the quark single particle energies. These findings have prompted us to think more 
thoroughly about the structure of baryons in such models and possible implications for 
the matter problem. We found that it takes only very little effort to generalize a previous 
Skyrme type treatment of the single, massless baryon to the case of baryonic matter. In 
the chiral limit, an extremely simple, yet non-trivial, picture emerges: Both the baryon 
and dense matter are described by a spatially varying chiral angle which is best charac- 
terized as a "chiral spiral" with constant helix angle. The number of windings within the 
full space of length L measures the number of baryons in the box. Since, by construction, 
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this kind of state does not cost any potential energy in addition to what is already stored 
in the vacuum, one does not have to pay the expected high price for delocalizing quarks. 
The energy density of baryonic matter is identical to that of a free Fermi gas of massless 
quarks, in spite of the presence of interaction effects. The same picture applies to the 
chiral Gross-Neveu as well as to the massless 't Hooft model and should be generic for 
all chiral large TV models. The baryon density is constant in space as a consequence of 
axial current conservation, and we have verified that the whole scenario is exact to lead- 
ing order in the 1/N expansion. In a slightly more speculative vein we then investigated 
modifications due to a small bare quark mass. Here our task was greatly facilitated by 
the fact that we only needed to pull together two independent investigations, the one of 



Ref. of the single baryon in field theoretic models with the one of Ref. |46[ of the 
sine-Gordon kink crystal, both inspired in some way by Skyrme's original ideas. As a re- 
sult, we have arrived at a rather comprehensive picture of matter at low and high density 
on the scale of the pion Compton wavelength. The crystal structure now becomes more 
conspicuous since also the baryon density displays a lattice of individual lumps. As an 
additional bonus, we have obtained a purely classical, mechanical model of what is going 
on (the sine-Gordon equation describes a chain of coupled pendulums, the quark mass 
playing the role of gravity). Given our starting point, namely the problem of baryonic 
matter in two dimensional quantum field theories like the Nambu-Jona-Lasinio model or 
QCD, this is rather amusing. 

It is noteworthy that a similar chiral structure of fermionic matter has been reported 
previously in a variety of models different from the present ones. This indicates that the 
basic results are more generally valid than our derivation might suggest. We mention here 
in particular the early work on the massive Schwinger model [0 and the more recent 



work on the massless Schwinger model with inert background charge |13| and QCD 2 with 
a finite number of colors and flavours jnj. Even more surprising are perhaps quite a 



number of speculations about spatially inhomogeneous chiral condensates with the same 
wave number as in our case, but in 3+1 dimensions, in the context of pion condensation 
18| , large iV QCD Hl5 , 16 , [17] , or effective chiral models [|l9j. In some of these works, 



the analogy with the Overhauser effect and spin-density waves (pairing of particle holes 
on opposite sides of the Fermi sphere) has been stressed. Although the language used is 
quite different from ours, there is no doubt that we are dealing with the same physical 
phenomenon. 

As a last remark, we wish to comment on the original Gross-Neveu model with 
only discrete chiral symmetry (pure (ipip) ^interaction). Most of the studies of the phase 
diagram for the Gross-Neveu model have in fact been performed for this model, and 
one might think that our analysis does not have anything to say about it. However, 
the criticism of Sect. 2 also applies here. Since the non-chiral Gross-Neveu model has 
(massive) bound baryons, the low density behaviour of the energy obtained in standard 
Hartree-Fock approximation cannot be correct, and the phase diagram may also have to 
be reconsidered. 



We thank A.C. Kalloniatis for a critical reading of the manuscript. Partial sup- 



16 



port of this work by the Bundesministerium fur Bildung, Wissenschaft, Forschung und 
Technologie is gratefully acknowledged. 



Appendix: Hartree-Fock single particle energies in the 
't Hooft model 

The Hamiltonian for the massless 't Hooft model in the axial gauge has the form 

H = Y, 2 -[(P+ V2) (a\(p) ai (p) - b}( P Mp)) + ^ E J;f(,/) ;¥" ,/) • (66) 

p,i ij,n^0 

Here we have regularized the theory by enclosing it in a box of length L with antiperiodic 
boundary conditions for the fermions. Note that this form is only valid in the limit L — > oo 
|22| . The CLi(p),bi(p) denote right- and left-handed quark operators, respectively, and the 



currents jij(n) can be taken in the U(iV) form at large N, 

jij( n ) = E ( a )(P) a i(P + n ) + b](p)bi(p + n)) . (67) 



It is important to understand that the Coulomb term still contains one- and two-body 
operators which can be disentangled by normal-ordering (up to 1/N corrections) as fol- 
lows, 

Y^kWd-n) = N J £(a\(p)a i (jp) + b\(p)b i (p)) (68) 

ij i,p 

- E ( a ](p) a j( ( l) a l( ( l + n ) a i(P + n ) + a ](p)bj(q)b\(q + n)ai(p + n) 



+b){p)aj(q)a\(q + n)bi(p + n) + &](p)&j(g)&J(g + n)bi(p + n) 

The Hamiltonian can be decomposed correspondingly into one- and two-body operators. 
Using the basic vacuum expectation values 

E<0|4(p)ai(?)|0) = f 8 m (l-cos9(p)) , 
£<0|&I(p)Mff)l°> = y^l+ooseCp)), 

i 

J2(0\a\(pMq)\0) = 52(0\bt(p)at(q)\0) = ~ 5 pq sm9(p) , (69) 
i i A 

the vacuum expectation value of these two contributions are found to be 



|#«|0> = NY,(-^(p+^cos9(p) + X<s2L 



, v L V 2J Kr/ 48 J 1 
(0\H^\0) = ivW-^£^[l + cos(%)-% + n))]] . (70) 

p \ JZ7r n ^0 U 
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The Ng 2 L/48 term in the 1-body part is just minus twice the "1" term in the 2-body 
part (remember that J2 n ^o = 7r 2 /3) . Hence, in the sum of both terms, the Coulomb 
energy involves the combination 

-i[cos(0(p)-0(p + n))-l] (71) 

where the infrared divergence has been tamed since denominator and numerator both 
vanish at n — 0. This cancellation between quark self-energy and Coulomb potential is 
similar to what happens in the meson equation of the 't Hooft model. The continuum limit 
then yields Eq. (EJ|). It is important to distinguish between one- and two-body operators 
here, because they enter with different relative weights in the single particle energies and 
in the total energy. Indeed, in the Hartree-Fock approach, if the single particle energies 
are decomposed according to their 1- and 2-body contributions as 

u(p) = J l \p) + uj {2 \ P ) , (72) 

then the ground state energy is 

(0\H\0)=Nj:^ 1) (p) + l^\p^ . (73) 

The factor 1/2 is necessary to avoid double counting of the 2-body interaction term. By 
comparison with Eq. ([7(]), we can turn this observation around and simply read off the 
single particle energies. We find in this way 

^ 2 \p) = -^e4( 1 + cos (^)-^ + ^))) • ( 74 ) 

iC)7r n^O 71 

Adding up the two contributions to ui(p), the "1" term is now cancelled instead of chang- 
ing sign. This is the reason why in the continuum limit we get the badly infrared divergent 
expression (^) for the quark energies. This short-cut derivation gives the same result as 
the more elaborate approach of Ref. [EH], where a single particle Hartree-Fock Hamil- 



tonian was first identified by commuting H with the quark operators and subsequently 
diagonalized. 
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Figure captions 

1. Physical fermion mass as a function of the Fermi momentum in the Gross-Neveu 
model, in units of m . 

2. Energy density per color as a function of the Fermi momentum in the Gross-Neveu 
model. Solid curve: Chirally symmetric solution (m = 0); diamonds: Broken chiral 
symmetry (m according to Fig. 1); dashed straight line: Mixed phase. Units of m . 

3. Quark condensate as a function of the Fermi momentum in the 't Hooft model, in 
units of (Ng 2 /2ny/ 2 . 

4. The complex condensate (ipip) +i(0i75 , 0) for the single baryon in units of (x/ji/j) v , as 
a function of x in units of L (chiral limit). 

5. Same as Fig. 4, but for baryonic matter. Each full turn of the spiral increases the 
baryon number by one unit. 

6. Same as Fig. 2 (Gross-Neveu model). Here we have included the energy density of 
the Skyrme crystal type of state (crosses), the true ground state. 

7. Solid curve: Spatial oscillation of the baryon density in the regime pf <C m n ; circles: 
Baryon density for a single baryon. 

8. Solid curve: Spatial oscillations of the scalar chiral condensate in the regime pj <C 
m n ; circles: Scalar chiral condensate for a single baryon. 

9. Same as Fig. 8, but for the pseudoscalar chiral condensate. 

10. Illustration of the distorted "chiral spiral" for baryonic matter at non-zero bare 
quark mass. 

11. Spatial dependence of baryon density as it evolves with increasing average density 
(or Fermi momentum), in units of m n . 
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